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Abstract 

We determine the 6-quark shape function parameters in the shape function scheme, mb(SF) and 
/i 2 (SF) using the Belle B — > X s ^ photon energy spectrum. We assume three models for the form 
of the shape function; exponential, gaussian and hyperbolic. 
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INTRODUCTION 



One of the best ways to determine the Cabbibo-Kobayashi-Maskawa(CKM) matrix ele- 
ment | V u b | is to measure the inclusive charmless semileptonic (B — > X u tv) partial branching 
fraction of the B meson. Its theoretical prediction relies on an accurate description of non- 
perturbative bound-state effects of a b quark in the B meson. The effects are encoded in 
.B-meson shape functions, which are not theoretically calculable and have to be determined 
experimentally. The leading order shape function can be determined using the photon en- 
ergy spectrum in B — > X s ^ decays, since up to the leading order of l/m& in the Heavy 
Quark Expansion it is described with the same leading order shape function as the one used 
for the prediction of the B — > X u iv decays. Such shape function determination was first 
performed by CLEO on their B — > X s ^ data [If, and was then repeated in Ref. for the 
Belle B — > X s ^ data, where the theory proposed by Kagan and Neubert [j| was used to 
describe the photon energy spectrum. 

Recently, Bosch, Lange, Neubert and Paz presented theoretical expressions for the triple 
differential B — > X u tv decay rates and for the B — > X s 7 photon spectrum, which incorporate 
all known contributions and smoothly interpolate between the "shape-function region" of 
large hadronic energy and small invariant mass, and the " OPE region" (in which all hadronic 
kinematical variables scale with Mb) [3, HI- The differential rate is given in terms of a leading 
shape function and a subleading shape function in the shape function scheme 0,1^, where in 
this scheme the leading shape function is expressed with the parameters mft(SF) and /i^(SF). 

In this paper we report on the determination of the shape function parameters by fitting 
the Belle B — ► X s j data with the predicted photon energy spectrum [5(, where the default 
subleading shape function model from Ref. p is used. It allows for a precise determination 
of \V u b\ as shown in Ref. and Ref. Q. 

Apart from using an updated and more complete theoretical description of the B — > 
X s 7 decays, this analysis differs from the one in Ref. by using a different set of shape 
function model parametrizations, namely exponential, gaussian and hyperbolic. It also 
includes subleading shape function effects that were absent in the former analysis. 



PROCEDURE 

We used a method based on that devised by the CLEO Collaboration We fit Monte 
Carlo (MC) simulated spectra to the raw data photon energy spectrum. "Raw" refers to 
the spectra that are obtained after the application of the B — > X s ^ analysis cuts. The use 
of "raw" spectra correctly accounts for Lorentz boost from the B rest frame to the center 
of mass system, energy resolution effects and avoids unfolding. The method is as follows: 

1. Assume a shape function model. 

2. Simulate the photon energy spectrum for a certain set of parameters; (m&(SF), 

/4(sf)). 

3. Perform a \ 2 fit of the simulated spectrum to the data where only the normalization 
of the simulated spectrum is floated and keep the resultant x 2 value. 

4. Repeat steps 2-3 for different sets of parameters to construct a two dimensional grid 
with each point having a \ 2 ■ 
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5. Find the minimum x 2 on the grid and all points on the grid that are one unit of \ 2 
above the minimum. 

6. Repeat steps 1-5 for a different shape function model. 



Shape function models 



Three shape function forms suggested in the literature are employed; exponential, gaus- 
sian and hyperbolic [5J]. Their functional forms are described in Table |U they are a function 
of Co and are parameterized by two parameters: A and b. Example shape function forms are 
plotted in Fig. 1. 
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TABLE I: The three models used for the shape function forms 



The parameters A and b are related to the HQET parameters A and ^ by analytical 
expressions Eq.46 and Eq.47 in Ref. [s| for exponential and gaussian models, respectively, 
while for the hyperbolic model the corresponding HQET parameters have to be calculated 
numerically. The shape function parameters mfc(SF) and /i^(SF) are obtained from the 
HQET parameters A and ^ using the relations in Eq.41 of Ref. (if, where the reference 
scale of 1.5 GeV is used. 



Monte Carlo simulated photon energy spectrum 

We generate B — > X s ^ MC events according to the prescription in Ref. for each set of 
the shape function parameter values. The generated events are simulated for the detector 
performance using the Belle detector simulation program and the B — > X s 7 analysis cuts 
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FIG. 1: Three different models with linear onset for small Cj. We show F( cxp ) (solid), F^ gauss ^ 
(dashed), and F^ hyp ) (dash-dotted), for parameters that correspond to m&(SF)= 4.63GeV/c 2 and 
/4(SF)= 0.2 GeV 2 /c 2 



are applied to the MC events to obtain the raw photon energy spectrum in the T(4S f ) rest 
frame 



Fitting the spectrum 

For a given set of shape function parameters, a x 2 fit of the MC simulated photon 
spectrum to the raw data spectrum is performed in the interval 1.8 < E * < 2.8 GeV [l^ ]. 
Although in the Ref. 0] the fitting was performed in the interval between 1.5 GeV and 
2.8 GeV, the data below 1.8 GeV are not used in the present analysis since the tails of the 
models we use below 1.8 GeV are not modelled accurately [Io| . 

The normalization parameter is floated in the fit. The raw spectrum is plotted in Figure 
the errors include both statistical and systematic errors. The latter are dominated by the 
estimation of the BB background and are 100% correlated. Therefore the covariance matrix 
is constructed as 

T r _stat „stat r i sys sys / , \ 

Vij = cr d . a d 6ij + a I a J (1) 
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where i,j = 1, 2, . . . , 10 denote the bin number, and a d is the error in the data. Then the 
X 2 used in the fitting is given by 
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where (V 1 )y denotes the ij th element of the inverted covariance matrix. The x 2 value of 
the fit is used to determine a map of x 2 as a function of the shape function parameters. 

The best fit and Ax 2 contour 



The best fit parameters are associated to the minimum chi-squared case, Xmm- The error 
"ellipse" is defined as the contour which satisfies A% 2 = x 2 ~ Xmm = 1- The contours are 
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FIG. 2: Raw B — > X s j photon energy spectra in the T(45) frame as acquired from data. The 
errors include both statistical and systematic errors. Raw refers to spectra as measured after the 
application of Belle B — > X s ^f analysis cuts. 



found to be well approximated by the function [llj, 

A^ (m ,(SF),,;(SF)) = ("'< SF ' + *<^' SF » a + + (<""< SF e » 2 + d ) 2 . (3) 

The parameters a, b, c, d, and e are determined by fitting the function to the parameter 
points that lie on the contour. 



RESULTS 



The best fit parameters are given in Table |HJ The parameter values are found to be 
consistent across all three shape function forms. The minimum x 2 fit for each shape function 
model is displayed in Figure El The fits to the contour with A% 2 = 1 points are shown in 
Figure 01 and |U The imposed shape function form acts to correlate mb(SF) and //^(SF). 
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TABLE II: The best ht shape function parameter values. 
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FIG. 3: Top: Minimum % 2 fits of MC simulated spectra to the raw data for each shape function 
model. Bottom: The fitted Ax 2 = 1 contours for each shape function model. 

SUMMARY 

The fe-quark leading shape function parameters in the shape function scheme, mb(SF) 
and /i^(SF), were determined from fits of Monte Carlo simulated spectra, generated by 
the prescription in Ref. |5j, to the raw Belle measured B — > X s 7 photon energy 
spectrum. Three models for the form of the leading shape function were used; expo- 
nential, gaussian and hyperbolic, while the default model from Ref. [j| was used for the 
subleading shape function, where the reference scale is chosen to be 1.5 GeV. Best fit 
parameters are: (m 6 (SF), /i 2 (SF)) cxp = (4.52,0.27), (m 6 (SF), /i 2 (SF)) gauss = (4.54,0.25), 
and (m fe (SF), /i 2 (SF)) hyp = (4.52,0.27), where m b (SF) and /z 2 (SF) are measured in units 
of GeV/c 2 and GeV 2 /c 2 respectively. We also determined the A% 2 = 1 contours in the 
(mf,(SF), /i 2 (SF)) parameter space for each of the assumed models. 
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FIG. 4: Comparison of the fitted Ax 2 = 1 contours for all shape function models. The contours 
for the exponential, gaussian and hyperbolic model are shown by the solid, dotted and dash-dotted 
curves, respectively. The vertical and horizontal lines mark the central values of the three fits. 
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